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A protocol is provided to reconstruct the Wigner function for the motional state of a trapped ion 
via fluorescence detection on another ion in the same trap. This "sympathetic tomography" of a 
dark ion without optical transitions suitable for state measurements is based on the mapping of its 
motional state onto one of the collective modes of the ion pair. The quantum state of this vibra- 
tional eigenmode is subsequently measured through sideband excitation of the bright ion. Physical 
processes to implement the desired state transfer and read-out are derived, and the accomplishment 
of the scheme for different mass ratios is evaluated. 
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I. INTRODUCTION 

In quantum physics, states are represented by a wave 
function or a density matrix, which may in turn be rep- 
resented as a phase space quasi-probability distribution. 
The Wigner function provides such a representation of 
the motional state of a particle, and it is convenient both 
because it provides a useful visualization of the position 
and momentum contents of the state [Ij and because it 
reveals non-classical properties of the system associated 
with quantum superposition states 

The ability to create and maintain superposition states 
is of relevance in quantum metrology and quantum in- 
formation applications, and it is a key element in at- 
tempts to investigate the quantum-classical correspon- 
dence through experimental studies of larger and more 
complex physical systems Experimental methods 

exist to reconstruct the motional state of a harmoni- 
cally trapped ion by subjecting it to laser fields and 
coupling its internal and motional degrees of freedom 
0, 13 ■ For particles with no known or useful optical 
transitions, however, reconstruction of the motional state 
would naturally proceed via monitoring of its spatial dis- 
tribution at different times after its release from the trap 
- a formidable task as it requires precise elimination of 
any fields that would perturb the ion motion and distort 
the imaging. 

In this paper we suggest to determine the motional 
quantum state of a trapped charged particle by use of 
an atomic ion in the same trap on which it is possible to 
perform laser excitation and fluorescence detection. We 
will refer to these particles as particle 1 and 2 and as the 
dark and bright ion in the following. The motion of the 
two ions is coupled via their Coulomb interaction, and 
our proposal is related to the idea of sympathetic cooling, 
where laser cooled ions exchange energy with other ions 
by their Coulomb interaction and thus bring them into 
thermal equilibrium at low temperature. 

The outline of the paper is as follows. In Sec. II, we 
determine the eigenmodes of the trapped ions and how 
they are related to the motion of the individual parti- 
cles. In Sec. Ill, we discuss how the Wigner function 




FIG. 1; (Color online) Schematic setup for sympathetic to- 
mography experiment. A dark and a bright charged particle 
are oscillating around their equilibrium positions in a har- 
monic one-dimensional ion trap. 



for the dark particle motion is related to the complete 
Wigner function for the two particle problem, and we 
show how a physical operation on the system makes the 
dark ion motion available by measurements of only one of 
the collective modes. In Sec. IV we analyze the physical 
parameters for different mass rations between the bright 
and the dark ion. In Sec. V we present a conclusion and 
outlook. 



II. TRAP EIGENMODES 

In this section, we investigate the axial vibrational 
modes of a pair of coupled ions in a linear Paul trap. 
This problem has been studied before, but unlike the 
conventional diagonalization of the potential in the La- 
grangian formulation, Q , we shall base our analysis on 
a diagonalization of the Hamiltonian into a pair of un- 
coupled oscillators. The eigenfrequencies and the precise 
form of the modes are equivalent in the Lagrangian and 
Hamiltonian approaches, but the concrete steps in our 
diagonalization of the Hamiltonian will become useful in 
the following. 

We assume that the ions are confined to a regime 
so that the trapping potential is well described by a 
quadratic expression in their excursions around their 
equilibrium position, and that they are transversally well 
confined, so that we can consider only the motion along 
the X-axis of a harmonic trap. 



The total Hamiltonian is then given by, 
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where rrii is the mass, and qi{pi) is the position (momen- 
tum) of the i*'* ion relative to its classical equilibrium 
, (« = 1, 2). In our work we con- 



position qi — X 
sider the situation where the dark ion is heavier than the 
bright ion, mi > m2. The internal state structure of ion 
1 is not specified, but we assume that ion 2 has a two- 
level internal structure that can be excited by resonant 
laser fields. 

In ([1]) , fc denotes the spring constant of harmonic mo- 
tion along the x-axis in the trap, related to the oscillation 
frequencies of a single ion of any of the two species by 
k = rriiUjf. 

The equilibrium distance between two ions of charge 

e is do = (iff^)^^^ 0' ^^'^ expanding the Coulomb 
interaction potential to second order the Hamiltonian of 
the system is given by. 
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where we use symbols qi and pi to recall that the posi- 
tions and momenta are operators, which fulfil the canon- 
ical communication relations, i.e. = ihSij. The 
Hamiltonian ^ has the form of two coupled harmonic 
oscillators, and it can conveniently be put in matrix form. 
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where, 
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We aim to find the transformation of the variables, 

X = MX, (6) 

leading to a diagonal Hamiltonian coefficient matrix 

B. = {A.r^fliM-\ (7) 

M must be a symplectic transform [l^ in order to 
preserve the canonical commutator relations. While a 
linear mixing of the two oscillators may bring the poten- 
tial terms represented in the upper left quadrant of H 
in diagonal form, in order to avoid the introduction of 
momentum cross terms in the lower right quadrant, we 
must first bring the kinetic energy operator into a diag- 
onal form. 



A. Squeezing transformation 

For this purpose, we introduce the scale, or squeezing, 
transformation, 
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where, s — (^)^^* is the rescaling parameter. This idea 
is illustrated and explained geometrically in [lo| . 

Under the transformation Eq. (|9]), the Hamiltonian 
coefficient matrix takes the form. 
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where /i = y/mim2. 



B. Mode mixing transformation 

At this point, one may perform a linear mixing of the 
two oscillator modes and eliminate the coupling terms 
in the potential. Note that this mixing is equivalent to 
a beam splitter transformation acting on two quantized 
fields, and it is equivalent to a spatial rotation of a two 
dimensional oscillator, conveniently parametrized by a 
rotation angle ip, 
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or in the compact form, 

X = M2X'. (12) 

With the correct choice of rotation angle, 
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it is possible to put the Hamiltonian into the desired 
diagonal form. 
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The eigenmode frequencies are given by, 




in agreement with the analysis in [9], and they represent 
collective motion of both ions. For the well studied case 
of identical ions, a;+ /cj_ = \/3, and the eigenmodes rep- 
resent the relative motion and center of mass motion of 
the ions. For different masses, the modes given by our ex- 
pressions are more complicated. The prospect of storing 
quantum states in the motional degree of freedom of one 
or several ions has led to numerous recent experiments 
which demonstrate the motional coupling in systems of 
identical and different ions The normal mode fre- 

quencies have thus been observed in experiments, both 
with trapped atomic ions such as Mg'^ and Be^ ^iB'], 
and as a very precise diagnostic tool, allowing the deter- 
mination of the mass of a dark molecular ion via probing 
of the collective mode sideband excitation frequency of a 
bright atomic ion 0, [l3| ■ 

The individual ion motion associated with the normal 
modes is identified by the final outcome of the squeezing 
Eq. ([5]) and the rotation transformations Eq. (fT^ . and 
is formally given by 

M = M2M1. (16) 

III. WIGNER FUNCTION RECONSTRUCTION 

In the previous section we showed that the normal 
modes of two different ions in a linear Paul trap are ob- 
tained by squeezing and rotation transformations of the 
single particle position and momentum coordinates. 

The complete Wigner function of both particles is a 
pseudo probability function W{qi,pi,q2,P2) on the four 
dimensional phase space of the particles. Alternatively, 
the Wigner function can be represented as a function 
of complex arguments W{ai, 02), where the real (imagi- 
nary) parts of ai are proportional to qi (pi). After convo- 
lution with a Gaussian, the Wigner function leads to the 
Husimi- or Q-function, which is the overlap of the state 
with the coherent (product) states \ai,a2)- 

In medical diagnostics, tomography uses a number 
of two-dimensional images recorded at different angles 
to reconstruct a three-dimensional image of a patient. 
Likewise, quantum state tomography reconstructs a two- 
dimensional phase space distribution by a mathematical 
transformation of sequentially measured one dimensional 
quadrature distributions. This method has been applied 
to traveling wave quantum fields [l^ , to the single mode 
field inside a cavity and to neutral and charged 

particle motion |24{ . 

For trapped ions, the reconstruction of the motional 
state is most conveniently carried out by a different 
method 0, which transfers information about the mo- 
tional quantum state onto the internal state of the ions, 
which is subsequently read out by fluorescence. If a 



trapped ion is excited with a laser on the lower side- 
band in the Lamb-Dicke regime, the Rabi frequency of 
excitation depends on the motional eigenstate via a cou- 
pling strength oc y^, and therefore monitoring of the 
internal state probabilities as function of time permits 
reconstruction of the motional eigenstate populations p„, 
Q. The |n = 0) vibrational ground state is equivalent to 
a zero amplitude coherent state, and its population po, 
thus represents the Q- function at the phase space origin. 
Displacement of the ion prior to detection permits detec- 
tion of the Q{a) function at any desired location a in 
phase space, and, moreover, as shown in Q, determining 
the full set of populations p„ for such displaced copies of 
a single quantum states leads, via a simple formula, di- 
rectly to the Wigner function W{a). For a single trapped 
ion, one applies the excitations on the lower motional 
sideband, and with two ions, we can use the difference in 
eigenfrequency to selectively address either of the eigen- 
modes and determine their separate Wigner functions 
W{a+) and W{a-). Note, however, that knowledge of 
W{a-\-) and W{a-) alone does not permit reconstruction 
of the full two-mode Wigner function, and also not of the 
Wigner function of the individual particles. 



A. Mapping on eigenmodes 

Whether one measures quadrature distribution func- 
tions or one assesses the phase space distribution via dis- 
placements and sideband excitations, the Wigner func- 
tion for two particles, two field modes or spatial motion 
in two dimensions, is a function of two complex or four 
real variables. The sampling of such a function is thus a 
very elaborate task. Fortunately, it is not needed to ac- 
quire information of the full two-mode Wigner function 
and extract from this quantity the required single parti- 
cle information. Instead, we propose to induce a physical 
process which maps the quantum state of particle 1 onto 
one of the eigenmodes, and subsequently perform single- 
mode tomography on that mode. 

What we seek is a physical process which maps the 
particle 1 state onto the (— ) oscillator state. This is 
equivalent to requesting that the unitary time evolution 
operator U (T) maps the particle 1 observables onto the 
similar observables for the (— ) mode. If the expansion of 
the Heisenberg picture operators q^{T),p^{T) coincides 
with the expansion of gi(0),pi(0) on the collective mode 
operators g±(0), _p±(0), measurements on the collective 
(— ) mode after the unitary mapping yields precisely the 
same outcome as one would obtain for the quadrature 
measurements on particle 1 in the initial state. The pre- 
vious section showed explicitly that the eigenmode linear 
combinations, indeed, are derived from the action of two 
physical operations: squeezing and mixing/rotation on 
the system, and we will now detail, which physical oper- 
ations, applied to the eigenmodes, will yield the desired 
mode mapping. 

The eigenmode operators and the individual particle 
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operators are related by the specific linear combinations 
expressed by the matrix M in Eg. dTBl) . with Mi and M2 
given in Eqs.(8,ll). Inverting, Eqs.(9,12), we can there- 
fore write the vector of individual particle quadrature 
operators in terms of the eigenmode operators 

X = (M2Mi)~^X" = M^^M^^X". (17) 

The derivation of the eigenmodes in Sec. II was not as- 
suming any actual implementation of the squeezing and 
mixing/rotation operations, but we note that by actu- 
ally performing such operations, the eigenmode operators 
evolve into, e.g., 

q-{T) = Ululq-{Q)U2Ui, (18) 

which is, indeed, given by the expression ([17]), if Ui 
squeezes <?" = q- and P2 = p+ by 1/s and their conju- 
gate variables by s (inverse of Eq.(8)), and U2 performs a 
rotation among the eigenmodes by the angle —ip (inverse 
of Eq.(ll)). 

In the Schrodinger picture, we wish to measure the 
time independent quadrature operators in the 

time evolved state of the system |*(r)) = U2Ui\^{Q)), 
or p{T) — U2Uip{0)UlU2, evolved first by the squeezing 
operation Ui on the eigenmodes and then by the rotation 
operation U2 among the eigenmodes of the system. 

B. Squeezing as a physical transformation 

In order to squeeze the collective modes of motion, we 
propose to follow the scheme in Ref. [6\ for motional 
squeezing of a single ion, by laser excitation detuned 
by twice the trap frequency. We address both collec- 
tive modes via the bright ion 2, and we use two classi- 
cal travelling waves with positive frequency parts with 
the same amplitudes and wave vectors n (or wave vector 
difference in case of Raman processes), but with differ- 
ent frequencies and phases, i?!^' = E^eA'^^^^'^'*''^'^ and 

E'j-p = £;oe'(''«2-""*-'^"). The plane wave fields refer 
directly to the position operator 52 of the illuminated ion. 
In order to squeeze the a = ± modes of collective mo- 
tion, the frequencies w/ = cueg + and w// = uJeg — 
are chosen, where Ueg is the frequency difference between 
the internal energy levels of the bright ion. Also, we put 
= —4>ii = — 7r/2. In the Lamb-Dicke regime, this 
effectively implements the Hamiltonian 

= zni],[77i2)]2(42 _ al)a^2, (19) 

(2) 

where the Lamb-Dicke parameters rja for the coupling 
of ion 2 to the (— ) and (+) modes depend on the nor- 
mal mode frequencies and on the masses of the two ions. 
Their values are provided in the Appendix, and they are 
shown in units of the Lamb-Dicke parameter Ty*-^-*, for a 
single trapped ion with mass TO2 in Fig. ([2]) as functions 
of the mass ratio toi/to2. As the mass ratio of the par- 
ticles increases, the Lamb-Dicke parameter for the cou- 

(2) 

pling of the bright ion to the high frequency mode 77Y 
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FIG. 2: (Color online) Lamb-Dicke parameters via for the 
coupling of the bright ion 2 to the eigenmodes. The Lamb- 
Dicke parameters are shown in units of the Lamb-Dicke pa- 
rameter ry'^' for a single ion with mass m2 as function of the 
mass ratio mi/m2. 

increases towards 0.8477*^^\ while its coupling to the lower 
frequency mode decreases towards zero. 

If the ion is prepared in the internal eigenstate {\g) + 
|e))/\/2 of the Pauli operator, Eq.(19) is a squeezing 
Hamiltonian for mode a, and we obtain the squeezing 
unitary operator of the mode a, 

f/i = e«"° (20) 

Ui provides the desired squeezing factor s when ^ = 
Q.,[7]'i^ft = ilns. Therefore, to squeeze the collective 
modes by the desired amounts, one has to adjust the in- 
teraction times such that, ta = hi s / {2Vla[ri^a'^Y) . The 
two upper (red and blue) curves in Fig. ([3]) show the 
dependence of the squeezing duration ta on the mass 
ratio in units of the single ion 2 squeezing time scale 
([?y(2)j2j-2-j-i_ ]\jQ^g th.aX for equal masses, no squeezing is 
required. 

C. Mixing/rotation as a physical operation 

After the squeezing has been applied to both mo- 
tional modes, we need to perform a rotation among the 
modes. This rotation is equivalent to a beam-splitter 
or mode- mixing transformation in optics 25], but the 
modes have different frequencies, so a simple constant 
coupling will not suffice. There exist, however, propos- 
als in cavity QED ^26,] to implement a mode mixing 
operation through the interaction of both modes with 
a single driven two-level system, and we shall adapt 
that scheme and illuminate the bright ion 2 with two 
classical standing wave lasers with positive frequency 
parts, Ef'^ = EoCOs{Kq2 + 0i)e-'^'^'*+'l"^ and E^f = 
EoCOs{Kq2 + 611)3-'^'^"*+'^"^ Here, both lasers have 
the same amplitudes Eo and wave numbers k, but dif- 
ferent frequencies and W77, which will enable the reso- 
nant transfer of excitation among the vibrational modes 
via the ion 2 internal state. We assume that the equi- 
librium position of the bright ion 2 is located at the 
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anti-node of both standing waves, i.e. 9i = On = 0. 
Under the above condition, the interaction Hamiltonian 
between the bright ion and standing wave lasers is given 

by, 

+ h.c], (21) 

where, Q. is the couphng Rabi frequency and a +2 = |e) {g\ 
is the raising operator of the internal ion state. The two 
laser frequencies are now chosen such that they bridge the 
frequency difference between the eigenmodes, a;/(//) = 

In the Lamb-Dickc regime, one may expand the in- 
teraction Hamiltonian of Eq. (piT) up to second order 
in ^2, and pass to the rotating frame given by the laser 
frequencies. Choosing 0/ — —cjiii — 7r/2, we obtain [26j . 

H = ihQ.i£\f^ {a^a\ - aUi+)cfx2, (22) 

Note that similar to our construction of the squeezing 
operators, we must prepare the ion in the internal super- 
position state \g) -f |e))/-\/2 to obtain the desired Hamil- 
tonian action on the oscillator modes, which in this case 
causes a perfect unitary mixing or rotation among the 
modes, 

[>(¥>)= e^(^-°+-^-°+\ (23) 

(2) (2) 

where, (p = flr]_ rj^ d:j;2t defines the rotation angle. To 
obtain the desired rotation angle f, given by Eq.(13), one 

(2) (2) 

has to adjust the interaction time such that ry^'ri_'VLt — 
Lf. The lower curve (green) in Fig. ^ shows this time as 
a function of the mass ratio mi/r7i2 in units of the single 
ion time scale ([?7'-^']^ri)~^, indicating that the duration 
of the rotation operation decreases with increasing mass 
ratio. 



D. Time scales 

In summary, we have presented physical operations 
that map the motional state of ion 1 onto the eigen- 
mode (— ), and we have shown that these operations can 
be carried out on time scales which are similar to the 
ones already demonstrated for the manipulation of single 
trapped ions. Since we selectively address the eigenmode 
sidebands, they must be resolved from the carrier and 
from each other, giving rise to the further requirement 
that the operation time must be sufficiently larger than 
(a;-|_ — a;_)~-^, w^^, and wZ^. In Fig. Q we show the 
interaction time for the mode mixing process in units 
of the inverse Rabi frequency (upper green curve) 
and compare it with (0;+ — a;_)~^, w^^, and wZ^ (lower 
curves). Note that the latter are shown in units of l/a;2, 
the inverse trapping frequency, and that long interaction 
times and resolved sidebands are guaranteed for both the 
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FIG. 3: (Color online) Duration of the mixing (green) and 
squeezing operations (red and blue) as functions of mxjm^. 
The durations are given in units of the single ion time scale 
([ry'^'j^r?)"^. Note that the duration of the squeezing of the 
high-frequency mode is multiplied by 10 for improved visibil- 
ity. 



squeezing and mixing processes if < ^2, which will 
quite generally be the case. 

In the beginning of this section we described how the 
Wigner function of the motional eigenmode can be de- 
termined by fluorescence measurements on the internal 
state of the bright ion, after displacement operations ap- 
plied to the eigenmode and excitation on the lower side- 
band transition. The diplacements may be done using 
the excitation mechanism provided in Ref. Q, and the 
fluorescence detection is efficiently carried out by means 
of the shelving method used also in Q . 

We note that the eigenmodes are not mixed by the 
Coulomb and trapping potential, and therefore the to- 
mographic measurement on the (— ) mode can be carried 
out in the interaction picture moving frame of that os- 
cillator with adequate time to perform measurements on 
the internal state of the bright ion. In practical applica- 
tions, however, the measurement time must of course be 
short compared to heating and decoherence time scales 
in the system. In the limit of heavy dark ions, the low 
frequency of the (— ) mode might make it more suscepti- 
ble to such decoherence, and one might instead transfer 
the state of particle 1 to the (-I-) mode with the higher 
frequency. This requires a larger rotation angle (/?, and 
hence the mapping process takes longer, suggesting that 
a criterion can be established for the optimum collective 
destination mode for the dark particle quantum state. 



IV. ANALYSIS 

In this section we present an overview of the quanti- 
tative estimates for our protocol in three different cases 
with Be^ as the bright ion. 
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FIG. 4: (Color online) Duration of mapping operations and 
fulfilment of sideband resolution. The upper (green) curve 
shows the duration of the mapping process in units of the 

(2) 

single ion time unit 1/Q, for = 0-2 vs. the mass ratio 
mi/m2. The lower curves show a;2/(i^+ — OJ-), UJ2/0J- and 
u)2/co+. As long as Q < ljJ2 the mapping is slow enough that 
the sidebands are well resolved. 
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FIG. 5: (Color online) Schematic setup for experiments with 
two ions trapped in separate local traps with harmonic oscilla- 
tor shape near the trapping minima. As the trap holding the 
bright ion is brought closer to the dark particle, the individual 
particle modes mix, and eigenmodes are formed. Using con- 
trolled displacement of the bright ion trap, squeezing of the 
eigenmodes, and mixing of the modes with bichromatic fields, 
the motional state of an initially isolated dark particle may 
be transferred to a collective mode and read out by sideband 
excitation. 



C. Larger mass differences 



A. Two ions with equal mass 

Consider the situation in which two Be'^ ions are lo- 
cated in the linear Paul trap. Making use of Eqs. (|Aip . 
the normal motional modes of the system are given by, 

q- = G.70gi +0.70^2, 

q+ = -0.7Ggi + 0.70^2, (24) 

with corresponding normal frequencies, a;_ = UJ2 and 
= \/3uj2 where W2 refers to the oscillation frequency 
of only ion 2 in the trap, lo- and a;+ correspond to 
the center-of-mass and stretch mode frequencies, and the 
normal modes are obtained by a mixing/rotation trans- 
formation with ip = 7r/4, c.f., Eq. In this case, the 

Lanib-Dicke parameters are evaluated as r;i^'' — 0.70r]^^^ 

and ry'^-' — 0.5477^^^ where 77*^^-' is the Lamb-Dicke param- 
eter of a single i?e+ ion located in the trap. 



B. A Be+ and an AI+ ion 

When a i?e+ ion is coupled to an Al^ ion, [1111/1112 = 3) 
and the normalized eigenmodes of motion are given by, 

q- = 0.98gi +0.21(j2, 

q+ = -0.54gi -1-0.84(72, (25) 

with corresponding eigenfrequencies, lo^ — 0.67w2 and 
a;+ = 1.49llI2. The low frequency in-phase mode of mo- 
tion has a larger share of the heavier particle's displace- 
ment, while particle 2 contributes a larger amount of the 
high frequency mode as witnessed by the Lamb-Dicke 
parameters, i}^^ = 0.4377^^) and ry!^^ = Q.77?7^^). 

In comparison to the case of equal masses, the mixing 
angle Lp = 0.36 is decreased and as a consequence the 
duration of the mixing process gets shorter. 



If the i?e+ ion, is coupled to a 20 times heavier particle 
in the trap, the collective normal modes correspond to, 

q- = 0.99gi + 0.02^2, 

q+ = -0.46gi-F 0.89^2, (26) 

with, w_ = 0.27a;2 and uj+ = 1.42^2, and ryi^^ = 0.2277(2) 
and vi"^^ ~ 0.837y(2'. The mixing angle p = 0.12, and 
by increasing the ratio 7711/7772 even further, the lower 
frequency converges to zero while the higher one tends 
to y/2uj2- 



V. CONCLUSION 

We have in this article proposed a method to recon- 
struct the Wigner function for the quantum state of mo- 
tion of a charged particle located in a linear ion trap, 
by transfer of the state to one of the eigenmodes shared 
with another ion in the trap on which laser excitation 
and fluoresecence detection is possible. The mapping of 
the quantum state is carried out by squeezing and mode 
mixing transformations which can all be implemented 
via illumination of the bright ion, and we showed that 
the operations are not more complex or time consum- 
ing than similar operations routinely applied to single 
trapped ions. 

The coupled motion of trapped ions has been stud- 
ied for a long time as a key element in applications for 
entanglement generation and quantum computing [2'^- 
[29| . where the motional state serves as a communication 
channel between the ion internal states. Recent investi- 
gations on the manipulation of the quantum states of mo- 
tion of several particles [T3l - [l5| show great potential for 
more delicate use of the motional degree of freedom, and 
we believe that the present work contributes to this new 
interesting research field. In particular, the demonstra- 
tion, following [26l| , that non-degenerate collective modes 
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can be effectively coupled by absorption and stimulated 
emission of fields with the appropriately chosen frequen- 
cies, may have much wider applications than indicated 
here. 

Our example dealt for simplicity with two ions located 
in the same trap, but with modern segmented surface 
electrode traps |30l], it is possible to prepare situations 
like the one depicted in Fig. 5, where two ions are located 
in separate harmonic potential minima. If the traps are 
far apart, the individual particle modes are not coupled, 
while for shorter distance, one can identify the eigen- 
modes and the associated squeezing and mixing trans- 
formations in the same manner as above. With time de- 
pendent potentials, it is also possible to move the traps, 
and bring the bright probe ion closer to the initially iso- 
lated dark particle. The oscillators with time dependent 
coupling are still governed by linear (symplectic) trans- 
formations of the quadrature observables and we can thus 
identify and carry out the physical transformations to 
map the quantum states for sympathetic tomography as 
we did in the single trap. The Bloch-Messiah reduction 
[Slj ensures that the desired transformation of the ob- 
servables can always be obtained by a combination of 
squeezing and mode mixing operations. 

Finally, we note the increasing interest in hybrid sys- 
tems, where, e.g., trapped ions may interact with neutral 
atoms [32:'| or molecules cold quantum gases 34 1, and 
optical cavity fields [s^ . The interactions are different in 
these cases, but similar methods may exist for sympa- 
thetic cooling and sympathetic tomography between the 
quantum states of the ion and of the partner system. 

Discussions with Michael Drewsen and Uffe V. Poulsen 
and financial support from the European Union Inte- 
grated Project AQUTE are gratefully acknowledged. 



applying rotation Eq. ^ and squeezing (ITT|) . we have. 



q- = s cos ipqi + s sm ipq2 , 
q+ = —ssmipqi+s^^cos(pq2, 



(Al) 



where, s — \Jm\lm2 and 93 is defined by Eq. (jl3p . In- 
verting the above equations, we obtain 



q2 — s sm ipq- 



scos ^pq+. 



(A2) 



Laser light couples to the internal degrees of freedom of 
ion 2 and thereby to the external degrees of freedom of 
the collective motion. In the dipole approximation, the 
coupling involves the plane wave factor e*''^^ , where n 
is the wave number of the light field (or the wave num- 
ber difference in case of a two-photon Raman process) 
and q2 is the position operator of the resonantly ex- 
cited particle. Substitution of q2 from Eq. (IA2p into 
the plane wave factor, and expansion of the collective 
mode position operators on rasing and lowering opera- 
nt ' 



tors, qa 



-{cLa 



leads to the indentification 



of the collective mode Lamb-Dicke factors. 



(2) . h 

V- = i^sirnp^ 



Appendix A 



Lamb-Dicke parameters 



(2) / ^ 



(A3) 



In order to identify the Lamb-Dicke parameters, one 
has to find the explicit form of the normal modes. After 



where the rotation angle is given by Eq.(13), and the 
mode eigenfrequencies are given by Eq.(15) in the text. 
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